Drosophila Models Uncover Substrate Channeling Effects on Phospholipids and Sphingolipids in Peroxisomal Biogenesis Disorders.

Peroxisomal Biogenesis Disorders Zellweger Spectrum (PBD-ZSD) disorders are a group of autosomal recessive defects in peroxisome formation that produce a multi-systemic disease presenting at birth or in childhood. Well documented clinical biomarkers such as elevated very long chain fatty acids (VLCFA) are key biochemical diagnostic findings in these conditions. Additional, secondary biochemical alterations such as elevated very long chain lysophosphatidylcholines are allowing newborn screening for peroxisomal disease. In addition, a more widespread impact on metabolism and lipids is increasingly being documented by metabolomic and lipidomic studies. Here we utilize Drosophila models of pex2 and pex16 as well as human plasma from individuals with PEX1 mutations. We identify phospholipid abnormalities in Drosophila larvae and brain characterized by differences in the quantities of phosphatidylcholine (PC) and phosphatidylethanolamines (PE) with long chain lengths and reduced levels of intermediate chain lengths. For diacylglycerol (DAG) the precursor of PE and PC through the Kennedy pathway, the intermediate chain lengths are increased suggesting an imbalance between DAGs and PE and PC that suggests the two acyl chain pools are not in equilibrium. Altered acyl chain lengths are also observed in PE ceramides in the fly models. Interestingly, plasma from human subjects exhibit phospholipid alterations similar to the fly model. Moreover, human plasma shows reduced levels of sphingomyelin with 18 and 22 carbon lengths but normal levels of C24. Our results suggest that peroxisomal biogenesis defects alter shuttling of the acyl chains of multiple phospholipid and ceramide lipid classes, whereas DAG species with intermediate fatty acids are more abundant. These data suggest an imbalance between de novo synthesis of PC and PE through the Kennedy pathway and remodeling of existing PC and PE through the Lands cycle. This imbalance is likely due to overabundance of very long and long acyl chains in PBD and a subsequent imbalance due to substrate channeling effects. Given the fundamental role of phospholipid and sphingolipids in nervous system functions, these observations suggest PBD-ZSD are diseases characterized by widespread cell membrane lipid abnormalities.

De novo variants in FRYL are associated with developmental delay, intellectual disability, and dysmorphic features.

FRY-like transcription coactivator (FRYL) belongs to a Furry protein family that is evolutionarily conserved from yeast to humans. The functions of FRYL in mammals are largely unknown, and variants in FRYL have not previously been associated with a Mendelian disease. Here, we report fourteen individuals with heterozygous variants in FRYL who present with developmental delay, intellectual disability, dysmorphic features, and other congenital anomalies in multiple systems. The variants are confirmed de novo in all individuals except one. Human genetic data suggest that FRYL is intolerant to loss of function (LoF). We find that the fly FRYL ortholog, furry (fry), is expressed in multiple tissues, including the central nervous system where it is present in neurons but not in glia. Homozygous fry LoF mutation is lethal at various developmental stages, and loss of fry in mutant clones causes defects in wings and compound eyes. We next modeled four out of the five missense variants found in affected individuals using fry knockin alleles. One variant behaves as a severe LoF variant, whereas two others behave as partial LoF variants. One variant does not cause any observable defect in flies, and the corresponding human variant is not confirmed to be de novo, suggesting that this is a variant of uncertain significance. In summary, our findings support that fry is required for proper development in flies and that the LoF variants in FRYL cause a dominant disorder with developmental and neurological symptoms due to haploinsufficiency.

Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma.

PURPOSE: YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases. METHODS: We report three unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs (GRs) of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6. RESULTS: Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type GR can rescue the lethality and autophagic flux defects whereas the variants are less efficient in rescuing the phenotypes. CONCLUSION: The YKT6 variants are partial loss-of-function alleles and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

Improving access to exome sequencing in a medically underserved population through the Texome Project.

PURPOSE: Genomic medicine can end diagnostic odysseys for patients with complex phenotypes; however, limitations in insurance coverage and other systemic barriers preclude individuals from accessing comprehensive genetics evaluation and testing. METHODS: The Texome Project is a 4-year study that reduces barriers to genomic testing for individuals from underserved and underrepresented populations. Participants with undiagnosed, rare diseases who have financial barriers to obtaining exome sequencing (ES) clinically are enrolled in the Texome Project. RESULTS: We highlight the Texome Project process and describe the outcomes of the first 60 ES results for study participants. Participants received a genetic evaluation, ES, and return of results at no cost. We summarize the psychosocial or medical implications of these genetic diagnoses. Thus far, ES provided molecular diagnoses for 18 out of 60 (30%) of Texome participants. Plus, in 11 out of 60 (18%) participants, a partial or probable diagnosis was identified. Overall, 5 participants had a change in medical management. CONCLUSION: To date, the Texome Project has recruited a racially, ethnically, and socioeconomically diverse cohort. The diagnostic rate and medical impact in this cohort support the need for expanded access to genetic testing and services. The Texome Project will continue reducing barriers to genomic care throughout the future study years.

Loss of the endoplasmic reticulum protein Tmem208 affects cell polarity, development, and viability.

Nascent proteins destined for the cell membrane and the secretory pathway are targeted to the endoplasmic reticulum (ER) either posttranslationally or cotranslationally. The signal-independent pathway, containing the protein TMEM208, is one of three pathways that facilitates the translocation of nascent proteins into the ER. The in vivo function of this protein is ill characterized in multicellular organisms. Here, we generated a CRISPR-induced null allele of the fruit fly ortholog CG8320/Tmem208 by replacing the gene with the Kozak-GAL4 sequence. We show that Tmem208 is broadly expressed in flies and that its loss causes lethality, although a few short-lived flies eclose. These animals exhibit wing and eye developmental defects consistent with impaired cell polarity and display mild ER stress. Tmem208 physically interacts with Frizzled (Fz), a planar cell polarity (PCP) receptor, and is required to maintain proper levels of Fz. Moreover, we identified a child with compound heterozygous variants in TMEM208 who presents with developmental delay, skeletal abnormalities, multiple hair whorls, cardiac, and neurological issues, symptoms that are associated with PCP defects in mice and humans. Additionally, fibroblasts of the proband display mild ER stress. Expression of the reference human TMEM208 in flies fully rescues the loss of Tmem208, and the two proband-specific variants fail to rescue, suggesting that they are loss-of-function alleles. In summary, our study uncovers a role of TMEM208 in development, shedding light on its significance in ER homeostasis and cell polarity.

Novel hemizygous single-nucleotide duplication in RPGR in a patient with retinal dystrophy and sensorineural hearing loss.

BACKGROUND: The RPGR gene has been associated with X-linked cone-rod dystrophy. This report describes a variant in RPGR detected with exome sequencing (ES). Genes like RPGR have not always been included in panel-based testing and thus genome-wide tests such as ES may be required for accurate diagnosis. METHODS: The Texome Project is studying the impact of ES in medically underserved patients who are in need of genomic testing to guide diagnosis and medical management. The hypothesis is that ES could uncover diagnoses not made by standard medical care. RESULTS: A 58-year-old male presented with retinitis pigmentosa, sensorineural hearing loss, and a family history of retinal diseases. A previous targeted gene panel for retinal disorders had not identified a molecular cause. ES through the Texome Project identified a novel, hemizygous variant in RPGR (NM_000328.3: c.1302dup, p.L435Sfs*18) that explained the ocular phenotype. CONCLUSIONS: Continued genetics evaluation can help to end diagnostic odysseys of patients. Careful consideration of genes represented when utilizing gene panels is crucial to ensure an accurate diagnosis. Medically underserved populations are less likely to receive comprehensive genetic testing in their diagnostic workup. Our report is an example of the medical impact of genomic medicine implementation.

De novo variants in PLCG1 are associated with hearing impairment, ocular pathology, and cardiac defects.

Phospholipase C isozymes (PLCs) hydrolyze phosphatidylinositol 4,5-bisphosphate into inositol 1,4,5-trisphosphate and diacylglycerol, important signaling molecules involved in many cellular processes. PLCG1 encodes the PLCγ1 isozyme that is broadly expressed. Hyperactive somatic mutations of PLCG1 are observed in multiple cancers, but only one germline variant has been reported. Here we describe three unrelated individuals with de novo heterozygous missense variants in PLCG1 (p.Asp1019Gly, p.His380Arg, and p.Asp1165Gly) who exhibit variable phenotypes including hearing loss, ocular pathology and cardiac septal defects. To model these variants in vivo, we generated the analogous variants in the Drosophila ortholog, small wing (sl). We created a null allele slT2A and assessed the expression pattern. sl is broadly expressed, including in wing discs, eye discs, and a subset of neurons and glia. Loss of sl causes wing size reductions, ectopic wing veins and supernumerary photoreceptors. We document that mutant flies exhibit a reduced lifespan and age-dependent locomotor defects. Expressing wild-type sl in slT2A mutant rescues the loss-of-function phenotypes whereas expressing the variants causes lethality. Ubiquitous overexpression of the variants also reduces viability, suggesting that the variants are toxic. Ectopic expression of an established hyperactive PLCG1 variant (p.Asp1165His) in the wing pouch causes severe wing phenotypes, resembling those observed with overexpression of the p.Asp1019Gly or p.Asp1165Gly variants, further arguing that these two are gain-of-function variants. However, the wing phenotypes associated with p.His380Arg overexpression are mild. Our data suggest that the PLCG1 de novo heterozygous missense variants are pathogenic and contribute to the features observed in the probands.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans.

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies.

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

Dicarboxylic acylcarnitine biomarkers in peroxisome biogenesis disorders.

The peroxisome is an essential eukaryotic organelle with diverse metabolic functions. Inherited peroxisomal disorders are associated with a wide spectrum of clinical outcomes and are broadly divided into two classes, those impacting peroxisome biogenesis (PBD) and those impacting specific peroxisomal factors. Prior studies have indicated a role for acylcarnitine testing in the diagnosis of some peroxisomal diseases through the detection of long chain dicarboxylic acylcarnitine abnormalities (C16-DC and C18-DC). However, there remains limited independent corroboration of these initial findings and acylcarnitine testing for peroxisomal diseases has not been widely adopted in clinical laboratories. To explore the utility of acylcarnitine testing in the diagnosis of peroxisomal disorders we applied a LC-MS/MS acylcarnitine method to study a heterogenous clinical sample set (n = 598) that included residual plasma specimens from nineteen patients with PBD caused by PEX1 or PEX6 deficiency, ranging in severity from lethal neonatal onset to mild late onset forms. Multiple dicarboxylic acylcarnitines were significantly elevated in PBD patients including medium to long chain (C8-DC to C18-DC) species as well as previously undescribed elevations of malonylcarnitine (C3-DC) and very long chain dicarboxylic acylcarnitines (C20-DC and C22-DC). The best performing plasma acylcarnitine biomarkers, C20-DC and C22-DC, were detected at elevated levels in 100% and 68% of PBD patients but were rarely elevated in patients that did not have a PBD. We extended our analysis to residual newborn screening blood spot cards and were able to detect dicarboxylic acylcarnitine abnormalities in a newborn with a PBD caused by PEX6 deficiency. Similar to prior studies, we failed to detect substantial dicarboxylic acylcarnitine abnormalities in blood spot cards from patients with x-linked adrenoleukodystrophy (x-ald) indicating that these biomarkers may have utility in quickly narrowing the differential diagnosis in patients with a positive newborn screen for x-ald. Overall, our study identifies widespread dicarboxylic acylcarnitine abnormalities in patients with PBD and highlights key acylcarnitine biomarkers for the detection of this class of inherited metabolic disease.

A comprehensive Drosophila resource to identify key functional interactions between SARS-CoV-2 factors and host proteins.

Development of effective therapies against SARS-CoV-2 infections relies on mechanistic knowledge of virus-host interface. Abundant physical interactions between viral and host proteins have been identified, but few have been functionally characterized. Harnessing the power of fly genetics, we develop a comprehensive Drosophila COVID-19 resource (DCR) consisting of publicly available strains for conditional tissue-specific expression of all SARS-CoV-2 encoded proteins, UAS-human cDNA transgenic lines encoding established host-viral interacting factors, and GAL4 insertion lines disrupting fly homologs of SARS-CoV-2 human interacting proteins. We demonstrate the utility of the DCR to functionally assess SARS-CoV-2 genes and candidate human binding partners. We show that NSP8 engages in strong genetic interactions with several human candidates, most prominently with the ATE1 arginyltransferase to induce actin arginylation and cytoskeletal disorganization, and that two ATE1 inhibitors can reverse NSP8 phenotypes. The DCR enables parallel global-scale functional analysis of SARS-CoV-2 components in a prime genetic model system.

A de novo missense variant in EZH1 associated with developmental delay exhibits functional deficits in Drosophila melanogaster.

EZH1, a polycomb repressive complex-2 component, is involved in a myriad of cellular processes. EZH1 represses transcription of downstream target genes through histone 3 lysine27 (H3K27) trimethylation (H3K27me3). Genetic variants in histone modifiers have been associated with developmental disorders, while EZH1 has not yet been linked to any human disease. However, the paralog EZH2 is associated with Weaver syndrome. Here we report a previously undiagnosed individual with a novel neurodevelopmental phenotype identified to have a de novo missense variant in EZH1 through exome sequencing. The individual presented in infancy with neurodevelopmental delay and hypotonia and was later noted to have proximal muscle weakness. The variant, p.A678G, is in the SET domain, known for its methyltransferase activity, and an analogous somatic or germline mutation in EZH2 has been reported in patients with B-cell lymphoma or Weaver syndrome, respectively. Human EZH1/2 are homologous to fly Enhancer of zeste (E(z)), an essential gene in Drosophila, and the affected residue (p.A678 in humans, p.A691 in flies) is conserved. To further study this variant, we obtained null alleles and generated transgenic flies expressing wildtype [E(z)WT] and the variant [E(z)A691G]. When expressed ubiquitously the variant rescues null-lethality similar to the wildtype. Overexpression of E(z)WT induces homeotic patterning defects but notably the E(z)A691G variant leads to dramatically stronger morphological phenotypes. We also note a dramatic loss of H3K27me2 and a corresponding increase in H3K27me3 in flies expressing E(z)A691G, suggesting this acts as a gain-of-function allele. In conclusion, here we present a novel EZH1 de novo variant associated with a neurodevelopmental disorder. Furthermore, we found that this variant has a functional impact in Drosophila.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder.

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

De novo variants in MRTFB have gain-of-function activity in Drosophila and are associated with a novel neurodevelopmental phenotype with dysmorphic features.

PURPOSE: Myocardin-related transcription factor B (MRTFB) is an important transcriptional regulator, which promotes the activity of an estimated 300 genes but is not known to underlie a Mendelian disorder. METHODS: Probands were identified through the efforts of the Undiagnosed Disease Network. Because the MRTFB protein is highly conserved between vertebrate and invertebrate model organisms, we generated a humanized Drosophila model expressing the human MRTFB protein in the same spatial and temporal pattern as the fly gene. Actin binding assays were used to validate the effect of the variants on MRTFB. RESULTS: Here, we report 2 pediatric probands with de novo variants in MRTFB (p.R104G and p.A91P) and mild dysmorphic features, intellectual disability, global developmental delays, speech apraxia, and impulse control issues. Expression of the variants within wing tissues of a fruit fly model resulted in changes in wing morphology. The MRTFBR104G and MRTFBA91P variants also display a decreased level of actin binding within critical RPEL domains, resulting in increased transcriptional activity and changes in the organization of the actin cytoskeleton. CONCLUSION: The MRTFBR104G and MRTFBA91P variants affect the regulation of the protein and underlie a novel neurodevelopmental disorder. Overall, our data suggest that these variants act as a gain of function.

Cluster Analysis of Short Sensory Profile Data Reveals Sensory-Based Subgroups in Autism Spectrum Disorder.

Autism spectrum disorder is a common, heterogeneous neurodevelopmental disorder lacking targeted treatments. Additional features include restricted, repetitive patterns of behaviors and differences in sensory processing. We hypothesized that detailed sensory features including modality specific hyper- and hypo-sensitivity could be used to identify clinically recognizable subgroups with unique underlying gene variants. Participants included 378 individuals with a clinical diagnosis of autism spectrum disorder who contributed Short Sensory Profile data assessing the frequency of sensory behaviors and whole genome sequencing results to the Autism Speaks' MSSNG database. Sensory phenotypes in this cohort were not randomly distributed with 10 patterns describing 43% (162/378) of participants. Cross comparison of two independent cluster analyses on sensory responses identified six distinct sensory-based subgroups. We then characterized subgroups by calculating the percent of patients in each subgroup who had variants with a Combined Annotation Dependent Depletion (CADD) score of 15 or greater in each of 24,896 genes. Each subgroup exhibited a unique pattern of genes with a high frequency of variants. These results support the use of sensory features to identify autism spectrum disorder subgroups with shared genetic variants.

De novo variants in FRMD5 are associated with developmental delay, intellectual disability, ataxia, and abnormalities of eye movement.

Proteins containing the FERM (four-point-one, ezrin, radixin, and moesin) domain link the plasma membrane with cytoskeletal structures at specific cellular locations and have been implicated in the localization of cell-membrane-associated proteins and/or phosphoinositides. FERM domain-containing protein 5 (FRMD5) localizes at cell adherens junctions and stabilizes cell-cell contacts. To date, variants in FRMD5 have not been associated with a Mendelian disease in OMIM. Here, we describe eight probands with rare heterozygous missense variants in FRMD5 who present with developmental delay, intellectual disability, ataxia, seizures, and abnormalities of eye movement. The variants are de novo in all for whom parental testing was available (six out of eight probands), and human genetic datasets suggest that FRMD5 is intolerant to loss of function (LoF). We found that the fly ortholog of FRMD5, CG5022 (dFrmd), is expressed in the larval and adult central nervous systems where it is present in neurons but not in glia. dFrmd LoF mutant flies are viable but are extremely sensitive to heat shock, which induces severe seizures. The mutants also exhibit defective responses to light. The human FRMD5 reference (Ref) cDNA rescues the fly dFrmd LoF phenotypes. In contrast, all the FRMD5 variants tested in this study (c.340T>C, c.1051A>G, c.1053C>G, c.1054T>C, c.1045A>C, and c.1637A>G) behave as partial LoF variants. In addition, our results indicate that two variants that were tested have dominant-negative effects. In summary, the evidence supports that the observed variants in FRMD5 cause neurological symptoms in humans.

Evidence for an association between Coffin-Siris syndrome and congenital diaphragmatic hernia.

Coffin-Siris syndrome (CSS) is an autosomal dominant neurodevelopmental syndrome that can present with a variety of structural birth defects. Pathogenic variants in 12 genes have been shown to cause CSS. Most of these genes encode proteins that are a part of the mammalian switch/sucrose non-fermentable (mSWI/SNF; BAF) complex. An association between genes that cause CSS and congenital diaphragmatic hernia (CDH) has been suggested based on case reports and the analysis of CSS and CDH cohorts. Here, we describe an unpublished individual with CSS and CDH, and we report additional clinical information on four published cases. Data from these individuals, and a review of the literature, provide evidence that deleterious variants in ARID1B, ARID1A, SMARCB1, SMARCA4, SMARCE1, ARID2, DPF2, and SMARCC2, which are associated with CSS types 1-8, respectively, are associated with the development of CDH. This suggests that additional genetic testing to identify a separate cause of CDH in an individual with CSS may be unwarranted, and that comprehensive genetic testing for individuals with non-isolated CDH should include an evaluation of CSS-related genes. These data also suggest that the mSWI/SNF (BAF) complex may play an important role in diaphragm development.

Complex effects on CaV2.1 channel gating caused by a CACNA1A variant associated with a severe neurodevelopmental disorder.

P/Q-type Ca2+ currents mediated by CaV2.1 channels are essential for active neurotransmitter release at neuromuscular junctions and many central synapses. Mutations in CACNA1A, the gene encoding the principal CaV2.1 α1A subunit, cause a broad spectrum of neurological disorders. Typically, gain-of-function (GOF) mutations are associated with migraine and epilepsy while loss-of-function (LOF) mutations are causative for episodic and congenital ataxias. However, a cluster of severe CaV2.1 channelopathies have overlapping presentations which suggests that channel dysfunction in these disorders cannot always be defined bimodally as GOF or LOF. In particular, the R1667P mutation causes focal seizures, generalized hypotonia, dysarthria, congenital ataxia and, in one case, cerebral edema leading ultimately to death. Here, we demonstrate that the R1667P mutation causes both channel GOF (hyperpolarizing voltage-dependence of activation, slowed deactivation) and LOF (slowed activation kinetics) when expressed heterologously in tsA-201 cells. We also observed a substantial reduction in Ca2+ current density in this heterologous system. These changes in channel gating and availability/expression manifested in diminished Ca2+ flux during action potential-like stimuli. However, the integrated Ca2+ fluxes were no different when normalized to tail current amplitude measured upon repolarization from the reversal potential. In summary, our findings indicate a complex functional effect of R1667P and support the idea that pathological missense mutations in CaV2.1 may not represent exclusively GOF or LOF.

Clinical diagnosis of metabolic disorders using untargeted metabolomic profiling and disease-specific networks learned from profiling data.

Untargeted metabolomics is a global molecular profiling technology that can be used to screen for inborn errors of metabolism (IEMs). Metabolite perturbations are evaluated based on current knowledge of specific metabolic pathway deficiencies, a manual diagnostic process that is qualitative, has limited scalability, and is not equipped to learn from accumulating clinical data. Our purpose was to improve upon manual diagnosis of IEMs in the clinic by developing novel computational methods for analyzing untargeted metabolomics data. We employed CTD, an automated computational diagnostic method that "connects the dots" between metabolite perturbations observed in individual metabolomics profiling data and modules identified in disease-specific metabolite co-perturbation networks learned from prior profiling data. We also extended CTD to calculate distances between any two individuals (CTDncd) and between an individual and a disease state (CTDdm), to provide additional network-quantified predictors for use in diagnosis. We show that across 539 plasma samples, CTD-based network-quantified measures can reproduce accurate diagnosis of 16 different IEMs, including adenylosuccinase deficiency, argininemia, argininosuccinic aciduria, aromatic L-amino acid decarboxylase deficiency, cerebral creatine deficiency syndrome type 2, citrullinemia, cobalamin biosynthesis defect, GABA-transaminase deficiency, glutaric acidemia type 1, maple syrup urine disease, methylmalonic aciduria, ornithine transcarbamylase deficiency, phenylketonuria, propionic acidemia, rhizomelic chondrodysplasia punctata, and the Zellweger spectrum disorders. Our approach can be used to supplement information from biochemical pathways and has the potential to significantly enhance the interpretation of variants of uncertain significance uncovered by exome sequencing. CTD, CTDdm, and CTDncd can serve as an essential toolset for biological interpretation of untargeted metabolomics data that overcomes limitations associated with manual diagnosis to assist diagnosticians in clinical decision-making. By automating and quantifying the interpretation of perturbation patterns, CTD can improve the speed and confidence by which clinical laboratory directors make diagnostic and treatment decisions, while automatically improving performance with new case data.